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Fusion of Multiple Imaging Planes for Isotropic Imaging in MRI and 
Quantitative Image Analysis Using Isotropic or Near-isotropic 
Imaging 



Statement as to Rights to Inventions Made Under Federally 
Sponsored Research and Development 

[0001] Certain aspects of the invention described below were made 
with United States Government support under Advanced Technology 
5 Program 70NANBOH301 6 awarded by the National Institute of Standards 
and Technology (NIST). The United States Government may have rights in 
certain of these inventions. 

Priority Claim 

[0002] This application claims the benefit of U.S. Provisional Patent 
10 Application Serial No. 60/431 ,176, filed December 4, 2002, the entirety of 
which is herein incorporated by reference. 

Technical Field 

[0003] This invention relates generally to medical imaging, and 

more specifically to medical imaging that facilitates analysis in more than 

15 one dimension, e.g. magnetic resonance imaging (MRI). More particularly 
the invention relates to isotropic imaging techniques used in medical 
imaging, such as MRI, to improve quantitative image analysis. 

Background of the Invention 
[0004] Magnetic resonance imaging (MRI) is a noninvasive imaging 

20 technique that provides clinicians and diagnosticians with information 

about the anatomical structure and condition of a region of interest within a 
subject. See, for example, U.S. Patent No. 5,671,741 to Lang et al. issued 
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September 30, 1997 for "Magnetic Resonance Imaging Technique for 
Tissue Characterization;" U.S. Patent No. 6,219,571 B1 to Hargreaves et 
al. issued April 17, 2002, for "Magnetic Resonance Imaging Using Driven 
Equilibrium Fourier Transform;" U.S. Patent No. 6,479,996 to Hoogeveen 
5 et al. issued November 12, 2002 for "Magnetic Resonance Imaging of 
Several Volumes;" U.S. Patent Application No. 2002/0087274 A1 to 
Alexander et al. published July 4, 2002 for "Assessing the Condition of a 
Joint and Preventing Damage." Commonly, in MRI, a substantially uniform 
temporally constant main magnetic field (B 0 ) is set up in an examination 

10 region in which a subject being imaged or examined is placed. Via radio 
frequency (RF) magnetic field (B-,) excitation and manipulations, selected 
magnetic dipoles in the subject that are otherwise aligned with the main 
magnetic field are tipped to excite magnetic resonance. The resonance is 
typically manipulated to induce detectable magnetic resonance echoes 

1 5 from a selected region of the subject. In imaging, the echoes are spatially 
encoded via magnetic gradients set up in the main magnetic field. The 
raw data from the MRI scanner is collected into a matrix, commonly known 
as k-space. By employing inverse Fourier, two-dimensional Fourier, three- 
dimensional Fourier, or other known transformations, an image 

20 representation of the subject is reconstructed from the k-space data. 

[0005] Conventional MRI scans produce a data volume, wherein the 

data volume is comprised of voxels having three-dimensional 
characteristics. The voxel dimensions are determined by the physical 
characteristics of the MRI machine as well as user settings. Thus, the 

25 image resolution of each voxel will be limited in at least one dimension, 
wherein the loss of resolution in at least one dimension may lead to three- 
dimensional imaging problems. 

[0006] There are many applications in which depth or three- 
dimensional ("3D") information is useful for diagnosis and formulation of 
30 treatment strategies. For example, in imaging blood vessels, cross- 
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sections merely show slices through vessels, making it difficult to diagnose 
stenosis or other abnormalities. Likewise, interventional imaging, such as 
needle tracking, catheter tracking, and the like, requires 3D information. 
Also, depth information is useful in the so-called interactive imaging 
5 techniques in which images are displayed in real or near-real time and in 
response to which the operator can adjust scanning parameters, such as 
view angle, contrast parameters, field of view, position, flip angle, 
repetition time, and resolution. 

[0007] Three-dimensional imaging generally involves either 

10 acquiring multiple two-dimensional or slice images that are combined to 
produce a volumetric image or, alternately, the use of three-dimensional 
imaging techniques. Much effort at improving the efficiency of volume 
imaging has been focused on speeding up the acquisition. For example, 
many two-dimensional fast scan procedures have been adapted to three- 
15 dimensional imaging. Likewise, efforts have been made to improve 
reconstruction speed and efficiency, for example, through the use of 
improved reconstruction algorithms. Nevertheless, three-dimensional 
imaging remains relatively slow. 

[00081 However, current MRI acquisition techniques do not provide 

20 high resolution in all planes and quantitative image analysis using isotropic 
or near-isotropic imaging. Accordingly, the present invention contemplates 
new and improved magnetic resonance imaging techniques. 

[0009] An additional problem not addressed by current 3D MRI 

scanning methods is the reduction of partial volume effects. Partial 
25 volume effects are caused when a voxel falls within the boundary between 
two scanned objects. For example, if a patient's knee is being sagittally 
scanned, a voxel may be orientated such that part of the voxel falls within 
the femur and part falls within a space outside of the femur. MR imaging 
will average the overall gray value over the entire voxel. The lower the 
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scanning resolution the greater the partial volume effects. In a 3D scan, 
where there is low resolution in at least one plane of the scan impact of the 
partial volume effects is greatly increased. Thus, there is a need for 
methods of forming 3D MRI scans with reduced impact of partial volume 
5 effects. 

[0010] Still further, an additional shortcoming of conventional 3D 

MRI scanning procedures is that boundaries of scanned objects may be 
missed due to scanning resolution and scan orientation. This may occur 
when a boundary of an object being scanned lies between the slice 
1 0 thickness of the scan or the boundary of an object is parallel to the 
imaging plane. Therefore there is a need for improved methods for 
reducing the likelihood of missed boundaries. 

Summary of the Invention 

[0011] The invention addresses the problem that with current 3D 
1 5 image acquisition techniques the in-plane (x-y plane) resolution of the 
slices is usually at least 3 times higher than the slice thickness (in z- 
dimension). The low resolution between the slices (typically in z-direction) 
leads to limitations with respect to 3D image analysis and visualization. 
The structure of 3-dimensional objects cannot be described with the same 
20 accuracy in all three dimensions. Partial volume effects affect 

interpretation and measurements in the z-dimension to a greater extent 
than in the x-y plane. Thus, resolution and accuracy of multiplanar 
reformations depend on the slicing direction through the volumetric data. 

[0012] In addition, the invention also addresses the issue of 

25 increasing accuracy of tissue segmentation and/or quantitative analysis of 
images, such as MR images. For example, after obtaining an isotropic or 
near-isotropic three-dimensional MR image {e.g., using pulse sequence 
acquisition techniques described herein and known in the field), particular 
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tissues can be extracted from the image with greater accuracy and, 
moreover are quantitative. Currently available subjective visual inspection 
techniques are not quantitative and, additionally, are often inaccurate. 

[0013] Thus, in one aspect, a method of improving resolution of 
5 images, such as MR images, is provided. In certain embodiments, the 
method includes, for example, obtaining at least two MR scans (e.g., 
scans in perpendicular planes) of a body part and merging the scans, 
thereby increasing resolution. In any of the methods described herein, the 
scans may be in any plane, for example, sagittal, coronal and/or axial 

1 0 imaging planes. Preferably, the second or subsequent scans contain a 
sufficient number of slices to cover the entire field of view of the first scan. 
Furthermore, in any of the methods described herein, the data obtained 
from the two or more scans are subsequently merged to form a new data 
volume, which is isotropic (or near-isotropic) and has a resolution 

1 5 corresponding to the in-plane resolution of S1 and S2. Merging may 
include, for example, determining a gray value for each voxel (V) of the 
new (merged) data volume. In certain embodiments, the gray values are 
obtained by: (a) determining the position in 3D space for V; (b) obtaining 
(e.g., from the original scans) gray values of the scans prior to fusion at 

20 this position; (c) interpolating (combining) gray values from S1 and S2 into 
a single gray value (G); and (d) assigning G to V. 

[0014] In any of the methods described herein, any living tissue can 

be imaged, including, but not limited to, joints, bones and/or organs (e.g., 
brain, liver, kidney, heart, blood vessels, GI tract, etc.). 
25 [0015] In accordance with the present invention there is provided a 

MRI scanning method, the method comprising, performing a first MRI scan 
of a body part in a first plane, wherein the first MRI scan generates a first 
image data volume; performing a second MRI scan of the body part in a 
second plane, wherein the second MRI scan generates a second image 



WO 2004/051301 



PCT/US2003/038682 



data volume; and combining the first and second image data volumes to 
form a resultant image data volume, wherein the resultant image data 
volume is isotropic. 

[0016] In accordance with another embodiment of the present 
5 invention there is provided a method for producing isotropic or near- 
isotropic image data, the method comprising: obtaining a first image data 
volume from a first MRI scan in a first plane; obtaining a second image 
data volume from a second MRI scan in a second plane; extracting 
boundary image data from each of the first and second image data 
10 volumes; combining said extracted boundary image data to form a 
resultant image data volume. 

[0017] In accordance with the present invention there is provided a 

method for generating a three dimensional data volume, the method 
comprising: acquiring at least two data volumes from at least two MRI 
1 5 scans performed in two different planes; combining the data volumes to 
form a resultant data volume; selecting a therapy in response to the 
resultant data volume; and deriving a shape for an implant. 

[0018] The system includes an image analysis method. The image 
analysis is performed by obtaining a first image of a body part in a first 

20 plane, wherein the first image generates a first image data volume; 

obtaining a second image of the body part in a second plane, wherein the 
second image generates a second image data volume; and combining the 
first and second image data volumes to form a resultant image data 
volume, wherein the resultant image data volume is isotropic. Additionally, 

25 first and second gray values can be obtained from the first and second 

image data volumes at one or more three-dimensional positions. That data 
can then be interpolated to provide a resultant gray value which is then 
assigned to a voxel in the three-dimensional position of the resultant data 
volume. As will be appreciated by those of skill in the art, the angle 
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between the images can range from about 0° to180°, or from 0° to 90°. 
Once these values have been obtained, a therapy or treatment can be 
selected to complement the data volume. A person of skill in the art will 
appreciate that at least one additional image of a body part taken in a 
5 plane different than any previous plane used can be taken used to 

generate additional image volume. From that image volume, data volume 
is generated which can then be combined with the first and second image 
data volumes to form a resultant data volume. Of course, extracting a 
boundary image data volume from the resulting image data volume can 
10 also be performed, if desired. 

[0019] A method is also described for producing isotropic or near- 
isotropic image data from images. This method generally comprises: 
obtaining a first image data volume from a first image in a first plane; 
obtaining a second image data volume from a second image in a second 

15 plane; extracting boundary image data from each of the first and second 
image data volumes; and combining the extracted boundary image data to 
form a resultant image data volume. Of course, it is possible to also 
obtaining at least one additional image data volume from at least one 
additional image in a plane different than the first plane and the second 

20 plane; extracting an additional boundary image data from the additional 
image data volume; and combining the additional boundary image data 
volume with the resultant image data volume. This resultant data can be 
isotropic or near-isotropic. As will be appreciated, the first plane can be at 
an angle relative to the second plane; that angle can be from about 0° to 

25 1 80° or from about 0° to 90°. 

[0020] A method is included for generating a three dimensional data 
volume. Generally, this method includes acquiring at least two data 
volumes from at least two images performed in two different planes; 
combining the data volumes to form a resultant data volume; selecting a 
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therapy in response to the resultant data volume; selecting an implant; and 
deriving a shape for an implant. The combining step can further include 
obtaining gray values for each data point in each of the data volumes; 
interpolating a resultant gray value from gray values; and assigning the 
5 resultant value to each data point of the resultant data volume. Prior to 
combining the data, data corresponding to any surface can be scanned in 
each plane and extracted. 

[0021] Another method for generating three dimensional data is also 
disclosed. This method includes obtaining a first image in a first plane 

10 producing a first data volume with a default resolution; obtaining a second 
image in a second plane producing a second data volume with the default 
resolution; combining the first and second data volumes to produce a 
resultant data volume, the resultant data volume having a resultant 
resolution. As will be appreciated, the resultant resolution is greater than 

15 the default resolution. 

[0022] An image analysis method is disclosed that includes the 

steps of obtaining at least one image of a body part in at least a first plane 
and a second plane, wherein the first plane generates a first image data 
volume and the second plane generates a second image data volume; and 

20 combining the first and second image data volumes to form a resultant 
image data volume, wherein the resultant image data volume is isotropic. 

[0023] An alternative image analysis method is also disclosed that 

includes obtaining at least one image of a body part in at least a first plane 
and a second plane, wherein the first plane generates a first image data 
25 volume and the second plane generates a second image data volume; and 
combining the first and second image data volumes to form a resultant 
image data volume, wherein the resultant image data volume is near- 
isotropic. 
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[0024] In accordance with the present invention there is provided a 

method for generating three dimensional MRI scan data, the method 
comprising: performing a first MRI scan in a first plane producing a first 
data volume with a default resolution; performing a second MRI scan in a 
5 second plane producing a second data volume with the default resolution; 
combining the first and second data volumes to produce a resultant data 
volume, the resultant data volume having a resultant resolution. 

Brief Description of the Figures 

[0025] Fig. 1 illustrates two MRI scans illustrating data volumes S1 
1 0 and S2; each of the scans shows a plurality of image slides taken in 
planes parallel to the initial scan. 

[0026] Fig. 2 illustrates a set of three voxels produced by an image 

scan illustrating an increased z-axis length. 

[0027] Fig. 3 illustrates a first set of three voxels produced by an 
1 5 image scan illustrating a z-axis component. 

[0028] Fig. 4 illustrates a second set of three voxels produced by an 
image scan illustrating a z-axis component. 

[0029] Fig. 5 illustrates a resultant set of nine voxels generated by 

the methods in accordance with the present invention. 

20 [0030] Fig. 6 illustrates a combined boundary image data extracted 
from two image scans. 

[0031] Fig. 7 illustrates a three-dimensional implant design 

generated from at least two image scans. 

[0032] Figs. 8a-c illustrate flow charts illustrating processes of the 
25 invention. 
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Detailed Description of the Preferred Embodiments 

[0033] The following description is presented to enable any person 
skilled in the art to make and use the invention. Various modifications to 
the embodiments described will be readily apparent to those skilled in the 
5 art, and the generic principles defined herein can be applied to other 

embodiments and applications without departing from the spirit and scope 
of the present invention as defined by the appended claims. Thus, the 
present invention is not intended to be limited to the embodiments shown, 
but is to be accorded the widest scope consistent with the principles and 
1 0 features disclosed herein. To the extent necessary to achieve a complete 
understanding of the invention disclosed, the specification and drawings of 
all issued patents, patent publications, and patent applications cited in this 
application are incorporated herein by reference. 

[0034] As will be appreciated by those of skill in the art, methods 
1 5 recited herein may be carried out in any order of the recited events which 
is logically possible, as well as the recited order of events. Furthermore, 
where a range of values is provided, it is understood that every intervening 
value, between the upper and lower limit of that range and any other 
stated or intervening value in that stated range is encompassed within the 
20 invention. Also, it is contemplated that any optional feature of the 

inventive variations described may be set forth and claimed independently, 
or in combination with any one or more of the features described herein. 
[0035] The present invention is of a method of image analysis that 

can be used for improving tissue segmentation andfor quantifying image 
25 analysis. Specifically, the present invention combines two or more images 
to achieve high resolution in all three-dimensional directions. The 
principles and operation of the method according to the present invention 
may be better understood with reference to the accompanying 
descriptions. 
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[0036] 1.0 General Overview 

[0037] According to the present invention, a method of improving 

resolution and/or tissue segmentation of images taken of a body part is 
described. This method typically involves acquiring at least two images in 
5 different planes and combining the images to achieve the same (e.g., high) 
degree of resolution in all directions. The images can be acquired, for 
example, by using an MRl. However, as other imaging devices become 
available, those of skill in the art will appreciate that these techniques can 
be provided to other imaging devices as well, without departing from the 

1 0 scope of the invention. 

[0038] The methods described herein provide isotropic or near- 

isotropic resolution which results in improved tissue segmentation. Unlike 
currently employed visual inspection, which is highly subjective, the 
methods and compositions described herein are quantitative and, 

1 5 accordingly, increase the accuracy of diagnosis and design of treatment 
regimes. 

[0039] 1.1 Magnetic Resonance Imaging (MRl) 

[0040] Describing MRl in general terms, all protons within living 

tissues have an inherent magnetic moment and spin randomly giving rise 

20 to no net magnetization or direction. When a specimen is placed within the 
magnetic field of the MR scanner, the protons continue to spin but align 
themselves parallel or anti-parallel to the direction of the field (B 0 ) 
corresponding to low and high-energy states respectively. In the course of 
an MR examination, a radiofrequency (RF) pulse (Bi) is applied to the 

25 sample from a transmitter coil orientated perpendicular to B 0 and the 
protons are momentarily tilted out of alignment; the precession of the 
induced net transverse magnetization around the axis of the static B 0 field 
produces a voltage across the ends of the receiver coil which is detected 
as the MR signal. For a general discussion of the basic MRl principles 
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and techniques, see MRI Basic Principles and Applications, Second 
Edition, Mark A. Brown and Richard C. Semelka, Wiley-Liss, Inc. (1999); 
see, also, U.S. Patent No. 6,219,571 to Hargreaves, et al. 

[0041] 1 .1 High Resolution 3D MRI Pulse Sequences 
5 [0042] MRI employs pulse sequences that allow for better contrast 

of different parts of the area being imaged. Different pulse sequences are 
better suited for visualization of different anatomic areas. More than one 
pulse sequence can be employed at the same time. A brief discussion of 
different types of pulse sequences is provided in International Patent 

1 0 Publication WO 02/22014 to Alexander et al. published March 21 , 2002. 
[0043] Routine MRI pulse sequences are available for imaging 

tissue, such as cartilage, include conventional T1 and T2-weighted spin- 
echo imaging, gradient recalled echo (GRE) imaging, magnetization 
transfer contrast (MTC) imaging, fast spin-echo (FSE) imaging, contrast 

15 enhanced imaging, rapid acquisition relaxation enhancement, (RARE) 
imaging, gradient echo acquisition in the steady state, (GRASS), and 
driven equilibrium Fourier transform (DEFT) imaging. As these imaging 
techniques are well known to one of skill in the art, e.g. someone having 
an advanced degree in imaging technology, each is discussed only 

20 generally hereinafter. 

[0044] 1.2. Measurement of T1 and T2 Relaxation 

[0045] As a result of random thermal motion, the proton spins within 

a sample lose coherence with one another. This loss of coherence results 
in signal decay. The time taken for the MR signal to return to zero depends 
25 on many factors, one is the rate at which the energized spins lose excess 
energy relative to their immediate environment. This phenomenon called 
spin-lattice, orT1 relaxation, affects mainly magnetization parallel to B 0 
and leads to a net loss of energy from the spin system. 
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[0046] Another phenomenon that is observed is that the spins of 
neighboring protons tend to drift out of alignment with one another as a 
result of slight differences in frequency. This causes a loss in phase 
coherence, referred to as spin-spin or T 2 relaxation. T 2 relaxation affects 
5 the transverse component of the magnetization but does not cause a net 
loss of energy. 

[0047] Conventional T1 and T2-weighted MRI depict living tissue 

such as articular cartilage, and can demonstrate defects and gross 
morphologic changes. One of skill in the art could readily select a T1 or 
10 T2-weighted MRI depending on the structure to be imaged. For example, 
T1 -weighted images show excellent intra-substance anatomic detail of 
certain tissue such as hyaline cartilage while T2-weighted imaging 
provides a better depiction of joint effusions and thus surface cartilage 
abnormalities. 

15 [0048] 1.3 Gradient-Recalled Echo (GRE) Imaging 

[0049] Gradient-recalled echo (GRE) imaging has 3D capability and 
the ability to provide high resolution images with relatively short scan 
times. Fat suppressed 3D spoiled gradient echo (FS-3D-SPGR) imaging 
has been shown to be more sensitive than standard MR imaging for the 
20 detection of hyaline cartilage defects such as those typically occurring in 
the knee. 

[0050] 1 .4 MAGNETIZATION TRANSFER CONTRAST IMAGING 

[0051] Magnetization transfer imaging can be used to separate 
articular cartilage from adjacent joint fluid and inflamed synovium. 
25 [0052] 1.5 Fast Spin-Echo (FSE) Imaging 

[0053] Fast spin-echo (FSE) imaging is another useful pulse 

sequence MRI technique. Incidental magnetization transfer contrast 
contributes to the signal characteristics of on fast spin-echo images and 
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can enhance the contrast between tissues. Sensitivity and specificity of 
fast spin-echo imaging have been reported to be 87% and 94% in a study 
with arthroscopic correlation. 

[0054] 1 .6 Echo Planar Imaging (EPI) 
5 [0055] Echo planar imaging (EPI) is an imaging technique in which 

a series of echoes is rapidly induced following a single radiofrequency 
(RF) pulse. More specifically, an RF pulse and a slice select gradient are 
applied to excite resonance in a selected slice and a phase encode 
gradient is applied to phase encode the resonance. A series of frequency 

10 encode or read gradients of alternating polarity is applied in successive 
fashion. During each read gradient, a magnetic resonance signal or echo 
is read out. Between each read gradient, a short pulse or blip along the 
phase encode gradient axis is applied to increment the phase encoding of 
the resonance by a line in the selected slice. A one-dimensional inverse 

1 5 Fourier transform of each echo provides a projection of the spin 

distribution along the read axis. A second inverse Fourier transform along 
the phase encoded echoes provides a second dimension of spatial 
encoding. Typically, the phase encode gradient blips are selected of an 
appropriate magnitude that data for a complete field of view is taken 

20 following each RF pulse. The total sampling time is determined by the 
number of sampled points per read gradient and the number of phase 
encode gradient steps. 

[0056] Echo volume imaging extends echo planar imaging 

techniques to multiple planes. After performing the above-described echo 
25 planar imaging sequence, a pulse or blip along a secondary phase 

encoding axis is applied. Typically, the secondary phase encoding blips 
step the phase encoding along an axis perpendicular to the primary phase 
encode and read axes. Thereafter, phase encode gradient blips are 
applied between each read gradient to step line by line in the primary 
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phase encode direction. Because the phase encode blips in the first k- 
space plane move the phase encoding to one extreme edge of the field of 
view, the phase encoding blips in the second k-space plane in the 
secondary phase encode direction are typically of the opposite polarity to 
5 step the phase encoding back in the opposite direction. In this manner, the 
multiple planes are aligned, but offset in steps in the z-direction. One 
disadvantage of the above echo planar imaging and echo volume imaging 
techniques is that the trajectory through k-space is reversed in time for 
alternate phase encode lines or views. This causes phase discontinuities 
10 that can result in ghosting. 

[0057] Spiral echo planar imaging techniques are also known, in 
which the applied x- and y-gradient pulses, i.e., along the traditional read 
and phase encode axes, are sinusoidally varying and linearly increasing. 
In this manner, data sampling commences at the center of the field of view 

1 5 and spirals outward, covering the field of view along a spiral k-space 

trajectory. One of the drawbacks of spiral echo planar imaging, however, is 
that it is a single slice technique. To obtain multiple slices, the spiral echo 
planar imaging technique is repeated multiple times. An RF excitation 
pulse and slice select gradient followed by sinusoidally varying and linearly 

20 increasing x and y-gradients are applied for each slice to achieve 
coverage of the volume of interest. 

[0058] 1.7 Contrast Enhancing Imaging 

[0059] The use of gadolinium in imaging has been applied in 

several different forms. For example, direct magnetic resonance (MR) 
25 arthrography, wherein a dilute solution containing gadolinium is injected 
directly into a tissue (e.g., joint), improves contrast between cartilage and 
the arthrographic fluid. Indirect MR arthrography, with a less invasive 
intravenous injection, can also been applied. Gadolinium enhanced 
imaging has the potential to monitor glycosaminoglycan content, which 
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may have implications for longitudinal evaluations of injured soft tissue 
such as cartilage. 

[0060] 1 .8 Driven Equilibrium Fourier Transformation 
[0061] Another 3D imaging method that has been developed is 
5 based on the driven equilibrium Fourier transform (DEFT) pulse sequence 
(U.S. Pat. No. 5,671,741 to Lang et al. issued September 30, 1997), and 
may be specifically utilized for soft tissue (e.g., cartilage) imaging. DEFT 
provides an effective tradeoff between T2/T1 weighting and spin density 
contrast that delineates the structures of interest. Contrast-to-noise ratio 

10 may, in certain tissues/structures, be greater with DEFT than with spoiled 
gradient echo (SPGR). DEFT is an alternative approach to SPGR. DEFT 
contrast is very well suited to imaging articular cartilage. Synovial fluid is 
high in signal intensity, and articular cartilage intermediate in signal 
intensity. Bone is dark, and lipids are suppressed using a fat saturation 

15 pulse. 

[0062] 1 .9 A Representative Example of MR Imaging 
[0063] A MR image can be performed using a whole body magnet 

operating at a field strength of 1.5 T (GE Signa, for example, equipped 
with the GE SR-120 high speed gradients [2.2 Gauss/cm in 184 usee 

20 risetimes]). Prior to MR imaging, external markers filled with Gd-DTPA 
(Magnevist.RTM., Berlex Inc., Wayne, N.J.) doped water (T1 relaxation 
time approximately 1 .0 sec) can be applied to the skin. External markers 
can be included in the field of view of all imaging studies. Patients can be 
placed in the scanner in supine position and the appropriate area imaged. 

25 After an axial scout sequence, coronal and sagittal T1 -weighted images 
can be acquired using the body coil (spin-echo, TR=500 msec, TE=15 
msec, 1 excitation (NEX), matrix 256 x128 elements, field of view (FOV) 
48 cm, slice thickness 7 mm, interslice spacing 1 mm). The scanner table 
can then be moved to obtain coronal and sagittal images using the same 
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sequence parameters. These T1 -weighted scans can be employed to 
identify axes that can be used later for defining the geometry of the tissue. 
A rapid scout scan can be acquired in the axial plane using a gradient 
echo sequence (GRASS, 2D Fourier Transform (2DFT), TR=50 msec, 
5 TE=1 0 msec, flip angle 40°, 1 excitation (NEX), matrix 256 x 1 28 
elements, field of view (FOV) 24 cm, slice thickness 7 mm, interslice 
spacing 3 mm). This scout scan can be used to determine all subsequent 
high resolution imaging sequences centered over the body part. 
Additionally, using the graphic, image based sequence prescription mode 
10 provided with the scanner software, the scout scan can help to ensure that 
all external markers are included in the field of view of the high resolution 
MR sequences. 

[0064] There are several issues to consider in obtaining a good 
image. One issue is good contrast between different tissues in the imaged 

1 5 area in order to facilitate the delineation and segmentation of the data 
sets. In addition, if there are external markers, these must be visualized. 
One way to address these issues is to use a three-dimensional spoiled 
gradient-echo sequence in the sagittal plane with the following parameters 
(SPGR, 3DFT, fat-saturated, TR=60 msec, TE=5 msec, flip angle 40°, 1 

20 excitation (NEX), matrix 256 x 160 elements, rectangular FOV 16 x 12 cm, 
slice thickness 1 .3 mm, 128 slices, acquisition time approximately 15 min). 
Using these parameters, one can obtain complete coverage across the 
body area and the external markers both in mediolateral and 
anteroposterior direction while achieving good spatial resolution and 

25 contrast-to-noise ratios. The fat-saturated 3D SPGR sequences can be 
used for rendering many tissues in three dimensions, e.g. cartilage. The 
3D SPGR sequence can then be repeated in the sagittal plane without fat 
saturation using the identical parameters and slice coordinates used 
during the previous acquisition with fat saturation. The resultant non-fat- 

30 saturated 3D SPGR images demonstrate good contrast between low 
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signal intensity cortical bone and high signal intensity bone marrow 
thereby facilitating 3D rendering of the femoral and tibial bone contours. It 
is to be understood that this approach is representative only for joints and 
should not be viewed as limiting in any way. 

5 [0065] 1.10 Magnetic Resonance Imaging -Vertically Open 
Magnet (0.5T) 

[0066] MR imaging can also be performed using a 0.5 T vertically 

open MR unit (GE Signa SP, General Electric, Milwaukee, Wis.) and a MR 
tracking system. Prior to MR imaging, external markers filled with Gd- 

10 DTPA (Magnevist.RTM., Berlex Inc., Wayne, N.J.) doped water (T1 
relaxation time approximately 1.0 sec) can be applied to the skin. The 
subject can be placed in upright position inside the magnet. The body part 
can be perpendicular to the main magnetic field. A 2DFT fast spin echo 
pulse sequence can be acquired in the sagittal plane (FSE, TR=4000 

1 5 msec, TE=25 msec, bandwidth 7.8 kHz, echo train length 8, 3 excitations, 
slice thickness 4 mm, interslice spacing 0.5 mm, matrix 256 x 192 
elements, field of view 24 cm). For rapid scan acquisition with scan plane 
tracking, a fast single slice gradient-echo pulse sequence can be acquired 
in the sagittal plane or in the axial plane (GRASS, TR=14 msec, TE=5 

20 msec, flip angle 40 degrees, bandwidth 32 kHz, 1 excitation, slice 

thickness 4 mm, matrix 256 x 128 elements, field of view 20 cm, temporal 
resolution 2 sec/image). A field of view of 20 cm can be chosen in order to 
achieve sufficient anatomic coverage in superoinferior. 
[0067] 2.0 Fusing Images 

25 [0068] Despite the existence of these imaging techniques, 

resolution in more than one plane remains difficult. In accordance with the 
present invention there is provided methods to overcome resolution 
difficulties wherein at least two data volumes from two separate images, 
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such as MRI scans, are combined to form a single data volume having 
isotropic or near-isotropic resolution. 

[0069] Referring now to Fig. 1, there is shown at least two 

exemplary data volumes S1 f 00 and S2 200 generated by at least two 
5 separate images. As illustrated here, each data volume 100, 200 has a 
plurality of data volumes 100 ( i. n} , 200 ( i. n) , as shown by the stacking of the 
image slices. S1 is an image of a knee joint taken in the coronal plane, 
while S2 is an image of a knee joint taken in the sagittal plane. In this 
example, S1 and S2 are taken in planes that are perpendicular to each 
10 other. However, as will be appreciated by those of skill in the art, other 
orientations and plane relationships can be used without departing from 
the scope of the invention. 

[0070] Each data volume can have equal imaging dimensions in two 

dimensions, for example, the x and y-axes, while the imaging dimension in 
15 a third dimension, e.g. z-axis, is greater than those in the first two 

dimensions, in this case the x-axis and y-axes. In a preferred embodiment 
a second scan can be taken at an angle between, for example, about 0° 
and 180° and more preferably between about 0° and 90°. 

[0071] Although the present invention is described using at least 

20 two scans, a person of skill in the art will appreciate that more scans can 
be used without departing from the scope of the invention. Thus, 
additional scans in the same or other planes or directions can also be 
obtained and analyzed. For example, if the first scan is acquired in the 
sagittal direction, a second scan in the coronal or axial imaging plane can 
25 then be acquired. 

[0072] It is possible that the second scan would have the same in- 

plane resolution as the first scan. The second scan could then contain a 
sufficient number of slices to cover the entire field of view of the first scan, 
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thereby resulting in two data volumes with information from the same 3D 
space. 

[0073] As described above, the data volumes generated from each 

image include x, y, and z-axis coordinate data, wherein as shown in Fig. 1 
5 , the x and y-axis data is isotropic while the z-axis data is non-isotropic. 
This phenomenon is better shown in Fig. 2 where there is shown an 
exemplary set of three voxels 200 as produced by an MRI scan in 
accordance with the present invention. The voxels 200 shown in Fig. 2 
are shown being orientated in the z-axis wherein the arrow 210 indicated 
10 the slice thickness of the image, in this case an MRI scan. The voxels 200 
further include a physical item 220 to be imaged. As shown in Fig. 2 it can 
be seen that due to the slice thickness 210, information pertaining to the 
physical item 220 to be imaged results in decreased accuracy. 

[0074] In addition to potentially missing data another problem with 
1 5 the larger slice thickness is the increase of partial volume effects. A 
partial volume effect occurs when a voxel only covers part of an object to 
be imaged, thus the gray value of the voxel is averaged instead of being a 
true gray value. As shown in Fig. 2, a partial volume effect occurs when a 
pixel or voxel is partially disposed over an object to be imaged 220. Since 
20 the voxel 200 is disposed partially over the object to be imaged the voxel's 
gray value will be averaged. To reduce the occurrence of partial volume 
effects, the present invention reduces the slice thickness of the scan, 
thereby reducing the likelihood of each voxel from being partially disposed 
on the object to be imaged. 

25 [0075] Referring now to Figs. 3-5, there is shown an exemplary 
embodiment of producing an isotropic or near-isotropic voxel in 
accordance with the present invention. As shown in Fig. 3, there is shown 
a set of three voxels 300 produced by an image scan, wherein the voxels 
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are shown having a z-axis 310 of greater length than the x-axis 315 and y- 
axis 320. 

[0076] Referring now to Fig. 4, there is shown a second set of three 

voxels 330 produced by a second image scan, wherein the second scan 
5 was taken at an angle 6 relative to the first scan. As discussed above, 0 
can range from, for example, about 0° to about 180°. As shown in Fig. 4, 
the second set of voxels has a z-axis dimension 340 greater than its x-axis 
345 and y-axis 350, wherein the z-axis 340 of the second set of voxels can 
be orientated in a plane different than that of the first set of voxels. 

10 [0077] Referring now to Fig. 5, there is shown a third set of voxels 
360 consisting of nine voxels, wherein the third set of voxels 360 has been 
formed by combining the first and second sets of voxels, wherein the z- 
axis data of the first set has been combined with x-axis or y-axis data from 
the second set of voxels to form a new z-axis 370 of the third set of voxels 

15 360, wherein the z-axis of the third voxel 370 has a length equal to or 
nearly equal to that of the x-axis 375 and y-axis 380, therefore producing 
voxels having isotropic or near-isotropic dimensions. 

[0078] After having performed at least two scans producing two 

data volume sets as shown in Fig. 1, the two data volumes are 

20 subsequently merged into a third data volume as shown in Figs. 3-5. This 
resultant data volume is isotropic or near-isotropic with a resolution 
corresponding to the in-plane resolution of S1 and S2. The gray value for 
each voxel of the third data volume is preferably calculated as follows: (a) 
determine the position in 3D space for each voxel; (b) determine (e.g., look 

25 up) the gray values in S1 and S2 at this position; (c) employ an 

appropriate interpolation scheme to combine the two gray values into a 
single gray value; and (d) assign each determined gray value to each 
voxel in the resultant data volume. 
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[0079] These manipulations may be repeated for more scans. 

Furthermore, to compensate for differences in positioning between Sf and 
S2 of the scanned subject, e.g. due to motion, a registration technique 
such as principal-axis or volume-based matching can be applied. 

5 [0080] 3.0 Tissue Segmentation 

[0081] In accordance with an alternative embodiment of the present 

invention there is provided a method of producing isotropic or near- 
isotropic MRI scan data from at least two image scans. 
[0082] As described in detail above, two individual data volumes are 

10 obtained from two separate image scans, wherein each of the scans have 
been taken at an angle 6 relative to each other. In a preferred 
embodiment, the second scan, or image, is preferably taken at an angle 6 
between about 0° and 180° more preferably between an angle between 
about 0° and 90°. Wherein each of the image scans produce individual 

15 data volumes having x, y, and z components wherein the x and y 
components are isotropic or nearly isotropic and the z-axis size is 
determined by the slice thickness (or step length) of, for example, the MRI 
machine. 

[0083] Tissue segmentation means can be applied to extract one or 
20 more tissues from one or more images. This can be achieved with 

classification of pixels or voxels of an electronic anatomical image (e.g. x- 
ray, CT, spiral CT, MRI) into distinct groups, where each group represents 
a tissue or anatomical structure or combination of tissues or anatomical 
structures or image background. For example, as described above, every 
25 data point of the first and second data volumes were combined to form a 
resultant data volume. Segmentation can then be performed on the entire 
data volume or subportions of the data volume. While effective in 
producing an isotropic or near-isotropic resultant data volume, the amount 
of data processing is great. Therefore, the method above requires a fairly 



WO 2004/051301 



PCT/US2003/038682 



23 

substantial amount of computer processing power as well as time to 
complete the mathematical calculations required. 

[0084] Referring now to Fig. 6 there is shown an exemplary 

embodiment of a three-dimensional MRI scan produced in accordance 
5 with the tissue segmentation methods of the present invention, wherein 
data pertaining to the object to be imaged has been first extracted from 
each of the data volumes prior to combining the extracted data to produce 
the three-dimensional scan of Fig. 6. As shown in Fig. 6, the cartilage 
surface of the medial femoral condyle is shown, wherein a sagittal scan 
10 400 and a coronal scan 450 were both acquired at a resolution of 0.27 mm 
x 0.27 mm in-plane resolution, and a 3mm slice thickness with a 0.5 mm 
spacing. As shown in Fig. 6, the medial edge is outlined well in the 
coronal scan, while the posterior edge of the condyle can be best seen in 
the sagittal scan. 

15 [0085] In accordance with the alternative embodiment, data 

pertaining to a surface or an area of interest is first extracted from each 
data volume produced by two or more image scans. After extracting the 
data volumes of interest, each pixel or voxel in the data volumes are 
subsequently merged into a new data volume by transforming them into a 

20 common coordinate system. This can be achieved by explicitly computing 
a transformation matrix for one or more of the data sets through the use of 
anatomical or other user defined landmarks or through a priori knowledge 
of the image position and orientation, such as the information provided by 
the DICOM imaging standard. To define a transformation in 3D space, the 

25 coordinates of four points in the original data volume and its corresponding 
location in the new data volume needs to be identified. These coordinate 
pairs are used to set-up a linear system of the form: 
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where A is the matrix with the original coordinates, T is the transformation 
matrix and B is the matrix with the new coordinates. The solution to the 
above system is given by: 

T = B x A" 1 

5 [0086] Alternatively, a transformation matrix can be implicitly 

calculated by performing a surface registration between the data sets. A 
surface registration algorithm merges the two data volumes by minimizing 
a cost function, such as a Euclidean distance transform and thus 
combining the volume data. Fig. 6 shows an example of a resultant data 
1 0 volume. This resultant data volume is isotropic or near-isotropic with a 
resolution corresponding to the in-plane resolution of S1 and S2. 

[00871 In another embodiment of the invention, a 3D MRI image is 

obtained using any suitable technique, for example using pulse sequence 
acquisition parameters that provide a 3D rather than a 2D Fourier 

1 5 Transform acquisition with isotropic or near-isotropic resolution, or by 
using fusion of two or more 2D acquisitions. As used herein, isotropic 
resolution refers to an MRI image in which the slice thickness is equal to 
the in-plane resolution. Similarly, the term "near-isotropic resolution" 
refers to an image in which the slice thickness does not exceed more than 

20 2x the in-plane resolution, more preferably not more than 1 .5x the in-plane 
resolution and even more preferably not more than 1.25x the in-plane 
resolution. The isotropic or near-isotropic 3DFT imaging pulse sequence 
has advantages with regard to partial volume averaging. Partial volume 
averaging is typically not greater in slice direction (z-direction) than in the 

25 imaging plane (x and y-direction). 

[0088] Non-limiting examples of pulse sequences suitable for 

obtaining near-isotropic or isotropic images include 3D FSE, 3D 
M FAST/3 D SS-SPGR, 3D FIESTA/3D SSFP, 3D FEMR, 3D DESS, 3D 
VIBE, and 3D SSFP. The preferred in-plane resolution of the 3DFT 
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isotropic or near-isotropic imaging sequence is less than 0.5mm and the 
preferred slice thickness is less than 0.8 mm, preferably less than 0.5 mm. 

[0089] Subsequently, these isotropic or near-isotropic resolution 

images are used to increase the accuracy of segmentation or tissue 
5 extraction and any subsequent visualizations and/or quantitative 

measurements of the body part, (e.g., measurement of cartilage thickness 
or size of cartilage defects). 

[0090] Thus, the invention described herein allows, among other 
things, for increased resolution and efficiency of tissue segmentation or 
10 tissue extraction. Following the manipulations described herein (e.g., 

merging of multiple images, isotropic or near-isotropic resolution imaging), 
commercially available segmentation software can be used, for example 
software that includes seed-growing algorithms and active-contour 
algorithms that are run on standard personal computers. 

15 [0091] For example, articular cartilage shown in the 3D MR images 
may be analyzed. A sharp interface is present between the high signal 
intensity bone marrow and the low signal intensity cortical bone thereby 
facilitating seed growing. 

[0092] One exemplary, but not limiting, approach uses a 3D surface 
20 detection technique that is based on a 2D edge detector (Wang-Binford) 
that has been extended to 3D. This surface detection technique can 
generate surface points and their corresponding surface normal. To 
smooth the contour, the program samples 25 percent of the surface points 
and fits a cubic spline to the sample points. The program can compute the 
25 curvature along sample spline points and find two sample points that have 
the maximum curvature and are separated by about half the number of 
voxels on the contour. These points partition the spline into two 
subcontours. For each subcontour, the program can compute the average 
distance between the points and the center of the mass. 
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[0093] Programs can allow the user, through the use of the mouse 
and/or keyboard, the ability to observe the scene from arbitrary angles; to 
start and stop the animation derived from the 3D data. Additionally, the 
user can derive quantitative information on the scene through selecting 
5 points with the mouse. 

[0094] The software programs can be written in the C++ computer 

language and can be compiled to run, for example on Silicon Graphics 
Workstations or Windows/Intel personal computers. 
[0095] 4.0 Three-dimensional Images 

1 0 [0096] After the 3D MR image is obtained, by either using a 3D 

acquisition or by fusing two or more 2D scans as described above, and 
after one or more anatomical objects have been extracted using 
segmentation techniques, for example, the object information can be 
transformed to a surface representation using a computer program. The 

15 program can, for example, be developed in AVS Express (Advanced 

Visual Systems, Inc., Waltham, Mass.). Every voxel has a value of zero if it 
is not within an object of interest or a value ranging from one to 4095, 
depending on the signal intensity as recorded by the 1 .5 T MR. An 
isosurface can then be calculated that corresponds to the boundary 

20 elements of the volume of interest. A tesselation of this isosurface is 

calculated, along with the outward pointing normal of each polygon of the 
tesselation. These polygons can be written to a file in a standard graphics 
format (e.g. Virtual Reality Modeling Language Version 1.0: VRML output 
language) and visualized on a computer screen. 

25 [0097] Visualization programs are also available, for example, user 
controllable 3D visual analysis tools. These programs read in a scene, 
which scene consists of the various 3D geometric representations or 
"actors." The program allows the user, through the use of the mouse 
and/or keyboard, the ability to observe the scene from arbitrary angles; to 
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start and stop the animation derived from the 3D data. Additionally, the 
user may derive quantitative information on the scene through selecting 
points with the mouse. 

[0098] The software programs can be written in the C++ computer 

5 language and be compiled to run, for example, on Silicon Graphics 
Workstations and Windows/Intel personal computers. Biochemical 
constituents, for example of cartilage, may also be visualized, for example 
as described in WO 02/22014 to Alexander. 

[00991 A method is also described for producing isotropic or near- 
1 0 isotropic image data from images. This method is shown in Fig. 8a. The 
first step is to obtain an image 800. As shown by optional repeat step 801, 
this step can be repeated such that multiple images are obtained. Suitable 
images include, for example, MRI. Once the image or images have been 
obtained image data volume is obtained 810. Image data volume can be 
1 5 obtained one or more times as indicated by the optional repeat step 811. . 
The generated image data is then combined to form an isotropic image 
volume 820 or a near-isotropic image volume 822. As will be appreciated 
by those of skill in the art, the process of forming one or more isotropic or 
near-isotropic image volume can be repeated one or more times as shown 
20 by optional repeat steps 821, 823. 

[0100] A method is also described for producing isotropic or near- 
isotropic image volume from images. This method is shown in Fig. 8b. 
The first step is to obtain an image data volume 830. As shown by optional 
repeat step 831, this step can be repeated such that multiple image data 
25 volumes are obtained. Once the image data volume or volumes have been 
obtained boundary image data is extracted 840„The extraction process 
can be repeated one or more times, as desired, 847. The extracted data 
volumes are combined to form isotropic image volumes 842, or near- 
isotropic image volumes 844. As will be appreciated by those of skill in the 
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art, more than one volume can be generated based on the extracted 
boundary image data 843, 845, if desired. 

[01 01] Another method is provided shown in Fig. 8c. This method 
includes the step of obtaining data volume from an image 850, a process 
5 that can optionally be repeated 851, if desired. Once the data volume is 
obtained, the data volume are combined to form at least one resultant data 
volume 860. More than one resultant data volume can be obtained 861, if 
desired. After obtaining the resultant data volume 860, a therapy can be 
selected based on the data volume 870 or a shape of an implant can be 
10 selected or derived 872. Either or both of these steps can be repeated, if 
desired, 871, 873. 

[01 02] Referring now to Fig. 7 there is shown an implant design 500 
generated by the methods in accordance with the present invention, 
wherein three-dimensional surface has been generated according the 

15 methods of the present invention. This three dimensional surface can then 
be utilized to manufacture an implant or select a therapy including an 
implant. Examples of such implants and implant techniques can be seen 

in co-pending US Patent Application No. __/__, filed on November 

25, 2003 by Aaron Berez, et al., for "Patient Selectable Joint Arthroplasty 

20 Devices and Surgical Tools Facilitating Increased Accuracy, Speed and 
Simplicity in Performing Total and Partial Joint Arthroplasty," the entirety of 
which is herein incorporated by reference. 

[0103] Alternatively, it may be determined after having generated a 
three-dimensional surface such as that shown in Fig. 7, that this type of 
25 implant is not necessary or cannot be utilized. Thus, a therapy may be 
chosen as an alternative or in addition to this type of implant. Examples of 
therapies include: drug therapy such as pain medication, 
chondroprotective agents, chondroregenerative agents, bone protective 
agents, bone regenerating agents, bone anabolic agents, bone osteoclast 
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inhibiting agents, injections of hyaluronic acid or chondroitin sulfates or 
other drugs or bioactive substances into the joint, osteotomy, chondral or 
osteochondral autografts or allografts, or other types of implants. Other 
types of implants can include, for example, total or unicompartmental 
5 arthroplasty devices. 

[0104] The methods present herein can be utilized at different 
timepoints, e.g. a timepoint T1 and a later or earlier timepoint T2. These 
timepoints can occur, for example, within one imaging session on a single 
day, or can occur over multiple imaging sessions over multiple days. The 
10 time span can further be hours, days, weeks, months and years. Tissues 
can then be characterized using quantitative measurements of the 
resultant data volumes V1 and V2 and changes in tissue composition or 
relative or absolute quantities can be assessed. 

[0105] The instant invention is shown and described herein in what 
15 is considered to be the most practical, and preferred embodiments. It is 
recognized, however, that departures may be made there from, which are 
within the scope of the invention, and that obvious modifications will occur 
to one skilled in the art upon reading this disclosure. 



WO 2004/051301 



PCT/US2003/038682 



30 
Claims 

WHAT IS CLAIMED IS: 

1 . An image analysis method comprising: 

obtaining a first image of a body part in a first plane, wherein 
the first image generates a first image data volume; 

obtaining a second image of the body part in a second plane, 
wherein the second image generates a second image data volume; 
and 

combining the first and second image data volumes to form a 
resultant image data volume, wherein the resultant image data 
volume is isotropic. 

2. The method according to claim 1 , wherein the combining 
step comprises: 

obtaining from said first and second image data volume first 
and second gray values at a three-dimensional position; 

interpolating a resultant gray value from said first and second 
gray values; and 

assigning said resultant value to a voxel at said three- 
dimensional position of said resultant data volume. 

3. The method according to claim 1 , wherein said second 
image is taken at an angle between about 0 and about 180 degrees from 
the first image. 

4. The method of claim 1 , wherein the first image is taken at a 
first angle and the second image is taken at a second angle, and further 
wherein the first angle does not equal the second angle. 

5. The method of claim 1 , wherein the first image is taken at a 
first time and the second image is taken at a second time. 
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6. The method of claim 3, wherein the second image is taken at 
an angle between about 0 and about 90 degrees from the first image. 

7. The method of claim 1 , further including: 

selecting a therapy in response to the resultant image data 

volume. 

8. The method of claim 1 , further including: 

selecting a treatment in response to the resultant image data 
volume. 

9. The method of claim 1 , further including: 

obtaining at least one additional image of a body part in a 
plane different than other planes, wherein the additional image 
generates an additional image data volume, wherein the additional 
data volume is combined with the first and second image data 
volumes to form a resultant data volume. 

1 0. The method of claim 1 , further including: 

extracting a boundary image data volume from the resulting image 
data volume. 

11. A method for producing isotropic or near-isotropic image 
data comprising: 

obtaining a first image data volume from a first image in a 
first plane; 

obtaining a second image data volume from a second image 
in a second plane; 

extracting boundary image data from each of the first and 
second image data volumes; and 
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combining the extracted boundary image data to form a 
resultant image data volume. 

12. The method of claim 1 1 , further including: 

obtaining at least one additional image data volume from at 
least one additional image in a plane different than the first plane 
and the second plane; 

extracting an additional boundary image data from the 
additional image data volume; and 

combining the additional boundary image data volume with 
the resultant image data volume. 

13. The method of claim 1 1 , wherein the resultant image data 
volume is near-isotropic. 

14. The method of claim 11, wherein the resultant image data 
volume is isotropic. 

15. The method of claim 1 1 , wherein the first plane is at an angle 
relative to the second plane. 

16. The method of claim 15, wherein the angle is between about 
0 and about 180 degrees. 

17. The method of claim 16, wherein the angle is between about 
0 and about 90 degrees. 

18. A method for generating a three dimensional data volume 
comprising: 

acquiring at least two data volumes from at least two images 
performed in two different planes; 

combining the data volumes to form a resultant data volume; 
selecting a therapy using the resultant data volume. 
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1 9. The method of claim 18, wherein the combining step 
comprises: 

obtaining gray values for each data point in each of the data 
volumes; 

interpolating a resultant gray value from gray values; and 
assigning the resultant value to each data point of the 
resultant data volume. 

20. The method of claim 19, further including; 

extracting data corresponding to a surface to be scanned in 
each of the planes before combining said data. 

21 . The method of claim 18, wherein the two scans are 
performed at ninety degrees relative to one another. 

22. A method for generating three dimensional data comprising: 
obtaining a first image in a first plane producing a first data 

volume with a default resolution; 

obtaining a second image in a second plane producing a 
second data volume with the default resolution; 

combining the first and second data volumes to produce a 
resultant data volume, the resultant data volume having a resultant 
resolution. 

23. The method of claim 22, wherein the resultant resolution is 
greater than the default resolution. 

24. An image analysis method comprising: 

obtaining at least one image of a body part in at least a first 
plane and a second plane, wherein the first plane generates a first 
image data volume and the second plane generates a second 
image data volume; and 
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combining the first and second image data volumes to form a 
resultant image data volume, wherein the resultant image data 
volume is isotropic. 

25. An image analysis method comprising: 

obtaining at least one image of a body part in at least a first 
plane and a second plane, wherein the first plane generates a first 
image data volume and the second plane generates a second 
image data volume; and 

combining the first and second image data volumes to form a 
resultant image data volume, wherein the resultant image data 
volume is near-isotropic. 

26. A method for generating a three dimensional data volume 
comprising: 

acquiring at least a first data volume and a second data 
volume from at least a first image and a second image, wherein the 
first image is obtained in a first plane and the second image is 
obtained in a second plane and further wherein the first plane is not 
equal to the second plane; 

combining the first data volume and the second data volume 
to form a resultant data volume; 

monitoring a therapy utilizing the resultant data volume. 

27. The method of claim 26, wherein the combining step 
comprises: 

obtaining gray values for each data point in each of the data 
volumes; 

interpolating a resultant gray value from gray values; and 
assigning the resultant value to each data point of the 
resultant data volume. 
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28. The method of claim 27, further including; 

extracting data corresponding to a surface to be scanned in 
each of the planes before combining said data. 

29. The method of claim 26, wherein the two scans are 
performed at ninety degrees relative to one another. 

30. The method of claim 26, wherein said treatment includes a 

drug. 

31 . The method of claim 26, wherein said treatment includes a 
bioactive agent. 

32. The method of claim 26, wherein said treatment includes 
surgery. 

33. The method of claim 26, wherein the first data volume is 
obtained at a first time point T1 and the second data volume is obtained at 
a second time point T2. 

34. A method for generating a three dimensional data volume 
comprising: 

acquiring at least two data volumes from at least two images 
performed in two different planes; 

combining the data volumes to form a resultant data volume; 
selecting an implant utilizing the resultant data volume. 

35. The method of claim 34, wherein the combining step 
comprises: 

obtaining gray values for each data point in each of the data 
volumes; 

interpolating a resultant gray value from gray values; and 
assigning the resultant value to each data point of the 
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resultant data volume. 

36. The method of claim 35, further including; 

extracting data corresponding to a surface to be scanned in 
each of the planes before combining said data. 

37. The method of claim 34, wherein the two scans are 
performed at ninety degrees relative to one another. 

38. A method for generating a three dimensional data volume 
comprising: 

acquiring at least two data volumes from at least two images 
performed in two different planes; 

combining the data volumes to form a resultant data volume; 
deriving an implant shape utilizing the resultant data volume. 

39. The method of claim 38, wherein the combining step 
comprises: 

obtaining gray values for each data point in each of the data 
volumes; 

interpolating a resultant gray value from gray values; and 
assigning the resultant value to each data point of the 
resultant data volume. 

40. The method of claim 39, further including; 

extracting data corresponding to a surface to be scanned in 
each of the planes before combining said data. 

41 . The method of claim 38, wherein the two scans are 
performed at ninety degrees relative to one another. 
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